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For the modelling of pedestrian dynamics we treat persons as self-driven objects moving in a 
continuous space. On the basis of a modified social force model we qualitatively analyze the influence 
of various approaches for the interaction between the pedestrians on the resulting velocity-density 
relation. To focus on the role of the required space and remote force we choose a one-dimensional 
model for this investigation. For those densities, where in two dimensions also passing is no longer 
possible and the mean value of the velocity depends primarily on the interaction, we obtain the 
following result: If the model increases the required space of a person with increasing current 
velocity, the reproduction of the typical form of the fundamental diagram is possible. Furthermore 
we demonstrate the influence of the remote force on the velocity-density relation. 

PACS numbers: 89.65.-s, 89.40.-a, 05.45.-a, 
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I. INTRODUCTION 



Microscopic models are state of the art for computer 
simulation of pedestrian dynamics. The modelling of the 
individual movement of pedestrians results in a descrip- 
tion of macroscopic pedestrian flow and allows e.g. the 
evaluation of escape routes, the design of pedestrian fa- 
cilities and the study of more theoretical questions. For 
a first overview see 0, The corresponding models 
can be classified in two categories: the cellular automata 
models p| , 0, IE an d models in a continuous space 
[cl I9L llOL lllj|. We focus on models continuous in space. 
They differ substantially with respect to the 'interaction' 
between the pedestrians and thus to the update algo- 
rithms as well. The social force model for example as- 
sumes, among other things, a repulsive force with remote 
action between the pedestrians jsl Ih l H^ . H I Hfl Hfil Il7j . 
Other models treat pedestrians by implementing a min- 
imum inter-person distance, which can be interpreted as 
the radius of a hard body jlOl [jjj . 

One primary test, whether the model is appropriate 
for a quantitative description of pedestrian flow, is the 
comparison with the empirical velocity-density relation 
[Til Il9l I20L l2lT | . In this context the fundamental dia- 
gram of Weidmann (2^ is frequently cited. It describes 
the velocity-density relation for the movement in a plane 
without bottlenecks, stairs or ramps. A multitude of 
causes can be considered which determine this depen- 
dency, for instance friction forces, the 'zipper' effect p3 | 
and marching in step pil |25| . As shown in [25J the em- 
pirical velocity-density relation for the single-file move- 
ment is similar to the relation for the movement in a 
plane in shape and magnitude. This surprising confor- 
mance indicates, that lateral interferences do not influ- 
ence the fundamental diagram at least up to a density- 
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value of 4.5 m~ 2 . This result suggests that it is sufficient 
to investigate the pedestrian flow of a one-dimensional 
system without loosing the essential macroscopic char- 
acteristics. We modify systematically the social force 
model to achieve a satisfying agreement with the em- 
pirical velocity-density relation (fundamental diagram). 
Furthermore we introduce different approaches for the 
interaction between the pedestrians to investigate the in- 
fluence of the required space and the remote action to 
the fundamental diagram. 



II. MODIFICATION OF THE SOCIAL FORCE 
MODEL 

A. Motivation 

The social force model was introduced by Q . It models 
the one-dimensional movement of a pedestrian i at posi- 
tion Xi (t) with velocity Vi (t) and mass m-i by the equation 
of motions 



dxi 
~~dt 



dvi - p - p ( 



(i) 



The summation over j accounts for the interaction with 
other pedestrians. We assume that friction at the bound- 
aries and random fluctuations can be neglected and thus 
the forces are reducible to a driving and a repulsive term 
Fj = Ff rv + F[ ep . According to the social force model 
we choose the driving term 
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drv 



mi 



(2) 



where v® is the intended speed and r» controls the ac- 
celeration. In the original model the introduction of the 
repulsive force F[ ep between the pedestrians is motivated 
by the observation that pedestrians stay away from each 
other by psychological reasons, e.g. to secure the private 
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sphere of each pedestrian 0. The complete model re- 
produces many self-organization phenomena like e.g. the 
formation of lanes in bi-directional streams and the os- 
cillations at bottlenecks @, O E3, E3, EH EH In 
the publications cited, the exact form of this repulsive 
interaction changes and the authors note that most phe- 
nomena are insensitive to its exact form [16j. We choose 
the force as in 13|. 



in order to meet the following properties: the force is al- 
ways pointing in the direction of the intended velocity v®; 
the movement of a pedestrian is only influenced by effects 
which are directly positioned in front; the required length 
d of a pedestrian to move with velocity v is d = a + bv. 

B. Interactions 



prep 



■VAidlxj-XiW-tk) 



-Bi 



(3) 



The hard core, di, reflects the size of the pedestrian i 
acting with a remote force on other pedestrians. Without 
other constraints a repulsive force which is symmetric in 
space can lead to velocities which are in opposite direc- 
tion to the intended speed. Furthermore, it is possible 
that the velocity of a pedestrian can exceed the intended 
speed through the impact of the forces of other pedes- 
trians. In a two-dimensional system this effect can be 
avoided through the introduction of additional forces like 
a lateral friction, together with an appropriate choice of 
the interaction parameters. In a one-dimensional sys- 
tem, where lateral interferences are excluded, a loophole 
is the direct limitation of the velocities to a certain in- 
terval HG3- 

Another important aspect in this context is the de- 
pendency between the current velocity and the space re- 
quirement. As suggested by Pauls in the extended ellipse 
model pij the area taken up by a pedestrian increase 
with increasing speed. Thompson also based his model 
on the assumption, that the velocity is a function of the 
inter-person distance 0. Furthermore Schreckenberg 
and Schadschneider observed in 0, , that in cellu- 
lar automata model's the consideration, that a pedes- 
trian occupies all cells passed in one time-step, has a 
large impact on the velocity-density relation. Hclbing 
and Molnar note in [j| that the range of the repulsive 
interaction is related to step-length. Following the above 
suggestion we specify the relation between required space 
and velocity for a one-dimensional system. In a one- 
dimensional system the required space changes to a re- 
quired length d. In [2jj it was shown that for the single- 
file movement the relation between the required lengths 
for one pedestrian to move with velocity v and v itself is 
linear at least for velocities 0.1 m/s < v < l.Qm/s. 



d = a+bv with a = 0.36 m and b— 1.06 s (4) 

Hence it is possible to determine one fundamental mi- 
croscopic parameter, d, of the interaction on the basis of 
empirical results. This allows focusing on the question 
if the interaction and the equation of motion result in a 
correct description of the individual movement of pedes- 
trians and the impact of the remote action. Summing up, 
for the modelling of regular motions of pedestrians we 
modify the reduced one-dimensional social force model 



To investigate the influence of the remote action both 
a force which treats pedestrians as simple hard bodies 
and a force according to Equation |3| where a remote ac- 
tion is present, will be introduced. For simplicity we set 
v® > 0, Xi+i > x% and the mass of a pedestrian to rrii = 1. 

Hard bodies without remote action 



W) = 



v°-Vi(t) 

-5(t)vi(t) 



x i+ i(t) - Xi(t) > di(t) 
x i+ i(t) - Xi(t) < di(t) 



(5) 



with 



di(t) = ai + biVi(t) 



The force which acts on pedestrian i depends only 
on the position, its velocity, and the position of the 
pedestrian i + 1 in front. As long as the distance between 
the pedestrians is larger than the required length, di, 
the movement of a pedestrian is only influenced by the 
driving term. If the required length at a given current 
velocity is larger than the distance the pedestrian stops 
(i. e. the velocity becomes zero). This ensures that 
the velocity of a pedestrian is restricted to the interval 
Vi = [0, v®] and that the movement is only influenced by 
the pedestrian in front. The definition of di is such that 
the required length increases with growing velocity. 

Hard bodies with remote action 



F(n -j : «i(*)>0 , fi x 

1 max (0,Gi(t)) : v % (t)<0 W 



with 



G((() = i^w_ e , 



1 



x i+1 (t) - Xi(t) - di(t) 



fi 



and 



diit) = a t + biVi(t) 



Again the force is only influenced by actions in front 
of the pedestrian. By means of the required length, di, 
the range of the interaction is a function of the velocity 
Vi. Two additional parameters, and fi, have to be 
introduced to fix the range and the strength of the 
force. Due to the remote action one has to change the 
condition for setting the velocity to zero. The above 
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definition assures that the pedestrian i stops if the force 
would lead to a negative velocity. With the proper 
choice of and /j and sufficiently small time steps 
this condition gets active mainly during the relaxation 
phase. Without remote action this becomes impor- 
tant. The pedestrian can proceed when the influence of 
the driving term is large enough to get positive velocities. 

This different formulation of the forces requires differ- 
ent update algorithms, which will be introduced in the 
next section. A special problem stems from the periodic 
boundary conditions enforced for the tests of the funda- 
mental diagram, as these destroy the ordering by causal- 
ity, which otherwise could avoid blocking situations. 

C. Time stepping algorithm 

The social force model gives a fairly large system of 
second order ordinary differential equations. For the hard 
body model with remote action, where the right hand 
side of the ODE's is continuous along the solution, an 
explicit Euler method with a time step of At = 0.001 s 
was tested and found sufficient. Within that time, the 
distance between two persons does not change enough to 
make the explicit scheme inaccurate. 

The situation for the hard body model without remote 
force is more complicated. Here the right hand side is a 
distribution, and the position of the Dirac spikes is not 
known a priory. Hence the perfect treatment is an adap- 
tive procedure, where each global time step is restricted 
to the interval up to the next contact. Unfortunately, 
this is a complicated and time consuming process. For 
a simple time step we choose the following procedure: 
Each person is advanced one step (At — 0.001 s) accord- 
ing to the local forces. If after this step the distance to 
the person in front is smaller than the required length, 
the velocity is set to zero and the position to the old po- 
sition. Additionally, the step of the next following person 
is reexamined. If it is still possible, the update is com- 
pleted. Otherwise, again the velocity is set to zero and 
the position is set to the old position, and so on. This is 
an approximation to the exact parallel update. It is not 
completely correct, however. To test its independence 
from the ordering of persons, computations using differ- 
ent orders were performed. The differences were minute 
and not more than expected from reordering of arith- 
metic operations. 



III. RESULTS 

To enable a comparison with the empirical fundamen- 
tal diagram of the single-file movement [25| we choose a 
system with periodic boundary conditions and a length 
of L = 17.3 m. For both interactions we proofed that 
for system-sizes of L = 17.3,20.0,50.0 m finite size ef- 
fects have no notable influence on the results. The values 



for the intended speed v® are distributed according to a 
normal-distribution with a mean value of /i = 1.24 m/s 
and a — 0.05 m/s. In a one-dimensional system the 
influence of the pedestrian with the smallest intended 
speed masks jamming effects which are not determined 
by individual properties. Thus we choose a er which 
is smaller than the empirical value and verified with 
er = 0.05, 0.1, 0.2 m/s, that a greater variation has no 
influence to the mean velocities at larger densities. 

In reality the parameters r, a, b, e and / are different 
for every pedestrian i and correlated with the individual 
intended speed. But we know from experiment ,25] that 
the movement of pedestrians is influenced by phenom- 
ena like marching in step and in reality the action of a 
pedestrian depends on the entire situation in front and 
not only on the distance to the next person. Therefore 
it's no point to attempt to give fully accurate values of 
this parameter and we may choose identical values for all 
pedestrians. We tested variations of the parameters and 
found that the behavior changes continuously. According 
to pjj , t = 0.61 s is a reliable value. 

For every run we set at t = all velocities to zero and 
distribute the persons randomly with a minimal distance 
of a in the system. After 3 x 10 5 relaxation-steps we 
perform 3 x 10 5 measurements-steps. At every step we 
determine the mean value of the velocity over all particles 
and calculate the mean value over time. The following 
figures present the dependency between mean velocity 
and density for different approaches to the interaction 
introduced in section Til Bl To demonstrate the influence 
of a required length dependent on velocity we choose dif- 
ferent values for the parameter b. With 6 = one get 
simple hard bodies. 
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FIG. 1: Velocity-density relation for hard bodies with a = 
0.36 m and without a remote action in comparison with em- 
pirical data from [25||. The filled squares result from sim- 
ple hard bodies. The introduction of a required length with 
b = 0.56 s leads to a good agreement with the empirical data. 
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Figure ^ shows the relation between the mean val- 
ues of walking speed and density for hard bodies with 
a = 0.36 m and without remote action, according to the 
interaction introduced in Equation [S] If the required 
length is independent of the velocity, one gets a nega- 
tive curvature of the function v — v{p). The velocity- 
dependence controls the curvature and b = 0.56 s re- 
sults in a good agreement with the empirical data. With 
b = 1.06 s we found a difference between the velocity- 
density relation predicted by the model and the empir- 
ical fundamental diagram. The reason for this discrep- 
ancy is that the interaction and equation of motion do 
not describe the individual movement of pedestrian cor- 
rectly. To illustrate the influence of the remote force, we 
fix the parameter a = 0.36 to, b — 0.56 s and set the val- 
ues which determine the remote force to e = 0.07 iV and 
/ = 2. 
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FIG. 2: Velocity-density relation for hard bodies with re- 
mote action in comparison with hard bodies without a re- 
mote action (filled circles). Again we choose a — 0.36m and 
b = 0.56 s. The parameter e = 0.07 N and f — 2 determine 
the remote force. With 6 = one gets a qualitative different 
fundamental diagram and a gap for the resulting velocities. 

The fundamental diagram for the interaction with re- 
mote action according to Equation [S] is presented in Fig- 
ure[2l The influence is small if one considers the velocity- 
dependence of the required length. But with 6 = one 
gets a qualitative different fundamental diagram. The in- 
crease of the velocity can be expected due to the effective 
reduction of the required length. The gap at p ~ 1.2 rn _1 
is surprising. It is generated through the development of 
distinct density waves, see Figure |31 as are well known 
from highways. From experimental view we have so far 
no hints to the development of strong density waves for 
pedestrians . The width of the gap can be changed by 
variation of the parameter / which controls the range of 
the remote force. Near the gap the occurrence of the den- 
sity waves depends on the distribution of the individual 
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FIG. 3: Time-development of the positions for densities near 
the velocity-gap, see Figure |2] For p > 1.2 m" 1 density waves 
are observable. Some individuals leave much larger than av- 
erage gaps in front. 



IV. DISCUSSION AND SUMMARY 

For the investigation of the influence of the required 
space and remote action on the fundamental diagram we 
have introduced a modified one-dimensional social force 
model. The modifications warrant that in the direction of 
intended speed negative velocities do not occur and that 
the motion of the pedestrians is influenced by objects and 
actions directly in front only. If one further takes into ac- 
count that the required length for moving with a certain 
velocity is a function of the current velocity the model- 
parameter can be adjusted to yield a good agreement 
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with the empirical fundamental diagram. This holds for 
hard bodies with and without remote action. The remote 
action has a sizeable influence on the resulting velocity- 
density relation only if the required length is independent 
of the velocity. In this case one observes distinct density 
waves, which lead to a velocity gap in the fundamental 
diagram. 

Thus we showed that the modified model is able to 
reproduce the empirical fundamental diagram of pedes- 
trian movement for a one-dimensional system, if it con- 
siders the velocity-dependence of the required length. 
For the model parameter b which correlates the required 
length with the current velocity, we have found that 
without remote action the value b — 0.56 s results in 
a velocity-density relation which is in a good agreement 
with the empirical fundamental diagram. However, from 
the same empirical fundamental diagram one determines 
b = 1.06 s, see [25|. We conclude that a model which re- 
produces the right macroscopic dependency between den- 
sity and velocity does not necessarily describe correctly 
the microscopic situation, and the space requirement of 
a person at average speed is much less than the average 
space requirement. This discrepancy may be explained 
by the 'short-sightedness' of the model. Actually, pedes- 



trians adapt their speed not only to the person imme- 
diately in front, but to the situation further ahead, too. 
This gives a much smoother movement than the model 
predicts. 

The above considerations refer to the simplest system 
in equilibrium and with periodic boundary conditions. 
In a real life scenario like a building evacuation, where 
one is interested in estimates of the time needed for the 
clearance of a building and the development of the den- 
sities in front of bottlenecks, one is confronted with open 
boundaries and conditions far from equilibrium. We as- 
sume that a consistency on a microscopic level needs to 
be achieved before one can accurately describe real life 
scenarios. The investigation presented provides a basis 
for a careful extension of the modified social force model 
and an upgrade to two dimensions including further in- 
teractions. 



Acknowledgments 

We thank Oliver Passon for careful reading and Wol- 
fram Klingsch for discussions. 



[1] M. Schreckenberg and S.D. Sharma (Ed.), Pedestrian and 

Evacuation Dynamics (Springer 2001) 
[2] E. R. Galea (Ed.), Pedestrian and Evacuation Dynamics 

(CMS Press 2003) 
[3] M. Muramatsu, T. Irie and T. Nagatani, Physica A267, 

487 (1999) 

[4] V. J. Blue and J. L. Adler, Journal of the Transportation 

Research Board 1678, 135 (2000) 
[5] K. Takimoto and T. Nagatani, Physica A 320, 611 

(2003) 

[6] C. Burstedde, K. Klauck, A. Schadschneider and J. Zit- 

tartz, Physica A 295, 507 (2001) 
[7] A. Kefiel, H. Kliipfel, J. Wahle and M. Schreckenberg, in 

Q, p. 193 

[8] D. Helbing and P. Molnar, Phys. Rev. E51, 4282 (1995) 
[9] S. P. Hoogendoorn and P. H. L. Bovy, in Proceedings 
of the 15th International Symposium on Transportation 
and Traffic Theory, Adelaide, 2002, edited by M. Taylor 
(University of South Australia, Adelaide, 2002), p. 219 
[10] P. Thompson and E. Marchant, Fire Safety Journal 24, 
131 (1995) 

[11] V. Schneider and R. Knnecke, in []J, p. 303 

[12] P. Molnar, Modellierung und Simulation der Dynamik 

von Fugangerstromen, Dissertation, (Shaker, Aachen, 

1996) 

[13] D. Helbing, I. Farkas and T. Vicsek, Phys. Rev. Lett. 84, 
1240 (2000) 



[14] D. Helbing, I. Farkas, and T. Vicsek, Nature 407, 487 
(2000) 

[15] D. Helbing, Rev. Mod. Phys. 73, 1067 (2001) 
[16] D. Helbing, I. Farkas, P. Molnar and T. Vicsek, in Q, 
p. 21 

[17] T. Werner and D. Helbing, in 0, p. 17 
[18] T. Meyer-Konig, H. Kliipfel and M. Schreckenberg, in Q, 
p. 297 

[19] A. Kirchner, H. Klpfel, K. Nishinari, A. Schadschneider 
and M. Schreckenberg, J. Stat. Mech. P10011 (2004) 

[20] S. P. Hoogendoorn, P. H. L. Bovy and W. Daamen, in 
Q, p. 123 

[21] www.rimea.de 

[22] U. Weidmann, Transporttechnik der Fufiganger, 

Schriftenreihe des IVT Nr. 90, zweite erganzte Au- 

flage, ETH Zurich (1993) 
[23] S. P. Hoogendoorn, W. Daamen, Transportation Science 

39/2, 0147 (2005) 
[24] P. D. Navin and R. J. Wheeler, Traffic Engineering 39, 

31 (1969) 

[25] A. Seyfried, B. Steffen, W. Klingsch and M. 
Boltes, available for download at http://arxiv.org 
physics/0506170 (2005) 

[26] J. L. Pauls, in Proceedings of the 3rd International Sym- 
posium on Human Behaviour in Fire, Belfast, 2004 (I n_ 
tersience Communications, London, 2004) 



